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Abstract 

We present a new approximation scheme for the centrifugal term to solve the Schrodinger equa- 
tion with the Hulthen potential for any arbitrary I state by means of a mathematical Nikiforov- 
Uvarov (NU) method. We obtain the bound state energy eigenvalues and the normalized corre- 
sponding eigenfunctions expressed in terms of the Jacobi polynomials or hyper geometric functions 
for a particle exposed to this potential field. Our numerical results of the energy eigenvalues are 
found to be in high agreement with those results obtained by using the program based on a nu- 
merical integration procedure. The s-wave (I = 0) analytic solution for the binding energies and 
eigenfunctions of a particle are also calculated. The physical meaning of the approximate analytical 
solution is discussed. The present approximation scheme is systematic and accurate. 
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I. INTRODUCTION 



It is necessary to obtain the exact bound state energy spectrum of the Schrodinger equa- 
tions for some physical potential models. Therefore, much works have been done to solve the 
wave equation for various radial and angular potentials. Unfortunately, the exact analytic 
solutions (EAS) of idealized quantum systems (QS), under consideration, are possible only 
in the s-wave case with angular quantum number I = for some exponential-type poten- 
tial models. On the other hand, the Schrodinger equation cannot be solved analytically 
for I 7^ because of the centrifugal term potential . Over the past years, some au- 
thors [1-16] have used the approximation l<yl ^ ~ ^~^y!2 for the centrifugal term potential 
proposed by Greene and Aldrich [1] to obtain the / ^ analytic bound-states [2,4,5] and 
scattering states [7] solutions of the non-relativistic [2,5] and relativistic [6] wave equations 
with some exponential-type potentials such as Hulthen potential [2-7], Eckart potential [10- 
13], Manning- Rosen potential [14-16] and diatomic molecular hyperbolical potential [17]. 
However, this approximation is valid only for small values of the screening parameter 5 and 
it breaks down for large values of 5 [5]. Therefore, there have been broad interest and im- 
pressive efforts in order to find a new approximation scheme which deals with the centrifugal 
term potential. 

The Hulthen potential [2,5,18] is the special case of the multiparameter exponential-type 
potential model [19,20]. It takes the form 

V(r) = -^ T , V = Ze*5, (1) 

where V is a constant and 5 is the screening parameter that determines the range of the 
potential. If the potential is used for atoms, then Vq = Z5 (in units h = c = e = 1), where 
Z is identified as the atomic number. The Hulthen potential behaves like the Coulomb 
potential near the origin (i.e., r — > 0) Vc(r) = —Ze 2 /r , but decreases exponentially in the 
asymptotic region when r ^> 0, so its capacity for bound states is smaller than the Coulomb 
potential [6,21-24]. This potential has been applied to a number of areas such as nuclear 
and particle physics [25-27], atomic physics [28-31], molecular physics [32,33] and chemical 
physics [34,35], etc. 

The bound-state EAS of the Schrodinger equation with the Hulthen potential can be 
solved in a closed form for s- waves (states with zero orbital angular momentum I) [36]. 
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However, for the case / 7^ 0, this quantum system cannot be exactly solved. For imple- 
menting approximate schemes economically and profitably; while dealing with practical 
quantum mechanical problems, EAS of the Hulthen potential is desirable although non- 
perturbative and numerical solutions of different potentials may lead to new physical ideas 
and/or calculational techniques in quantum physics. For instance, the numerical integration 
of the Schrodinger equation [37] is used to obtain the energy eigenvalues numerically for the 
Hulthen potential case, this provides a probe and/or test for the exactness of any analytic 
solution. One-parameter variational calculations are carried out in such numerical integra- 
tion methods. The variational results are practically identical to the exact energies, except 
in the high-screening region. These variational calculations turn to become sophisticated in 
the solution of Schrodinger equation with multi-parameter potentials. However, no "exact" 
values obtained from a numerical integration of the Schrodinger equation have been available 
to assess the accuracy of the various methods [37]. Hence, it is important to note that the 
analytic solution of any quantum potential model, even if it is an approximated solution, 
is indispensable since the obtained expressions for energy eigenvalues and eigen functions 
contain all the necessary information regarding the quantum system under consideration. 

In the non-relativistic case, for / 7^ 0, several techniques have been used to obtain approx- 
imate analytic solutions, some authors have obtained the bound-state energy eigenvalues by 
using the numerical integration approach [37,38], quasi-analytical variational [37,39], pertur- 
bation [40], SUSYQM [3], shifted 1/N expansion [41], AIM [5] and Nikiforov-Uvarov (NU) 
[2] methods. The results obtained by some of these methods [3,5] are in good agreement 
with the numerical integration approach [37] for low-screening region (small values of the 
screening parameter 5) but the agreement becomes poor in the high-screening regime [5]. 

Recently, Haouat and Chetouani [42] have solved the Klein-Gordon and Dirac equations 
in the presence of the Hulthen potential, where the energy spectrum and the scattering 
wavefunctions are obtained for spin-0 and spin-| particles, making a slight modification 
to the usual approximation scheme, ^ ~ a 2 ^ e _^ r ^ where 7 is a dimensionless parameter 
(7 = 0,1 and 2) for the centrifugal term potential. They found that the good approximation, 
however, when the screening parameter a and the dimensionless parameter 7 are taken as 
a — 0.1 and 7 =1, respectively, which is simply the case of the normal approximation used 
in the literature. Also, Jia and collaborators [43] have recently proposed a new alternative 




where a; is a dimensionless parameter 
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(cu = 1.030), for the centrifugal potential to improve the numerical energy eigenvalues of the 
Hulthen potential. When taking uj — 1, their approximation can be reduced to the usual 
approximation [1-16]. However, the accuracy of their numerical results [43] is still in poor 
agreement with the other numerical integration and variational methods [37] especially in 
high-screening S regime. This problem could be solved by making a better approximation 
for the centrifugal term potential. In this work, for any arbitrary /-state, we aim to ob- 
tain approximate energy eigenvalues and corresponding normalized wave functions for the 
Hulthen potential in high agreement with the numerical method [37]. Hence, we present an 
alternative effective approximation that gives highly accurate numerical energy eigenvalues 
of the Hulthen potential as a function of screening parameter for all states with Z — 1. 

This paper is organized as follows: In the next Section, the NU method is briefly intro- 
duced. In Section 3, the /-states Schrodinger equation for the Hulthen potential is solved 
within the new effective approximation scheme and using the NU method. The calculated 
energy eigenvalues and wave functions are compared with the other ones found by using 
different analytical and numerical methods. The normalized wave functions are obtained in 
Section 4. Finally, the relevant conclusions are given in Section 5. 



II. NU METHOD 



The NU method is briefly introduced here and the details can be found in Nikiforov- 
Uvarov handbook [44]. This method was proposed to solve the second-order differential 
wave equation of the hypergeometric-type: 

+ ^/n(s) + ^M*) = 0, (2) 

where <r(s) and a(s) are polynomials, at most of second- degree, and r(s) is a first-degree 
polynomial. The prime sign denotes the differentiation with respect to s. To find a particular 
solution of Eq. (2), one can decompose the wavefunction ip n (z) as follows: 

^n(s) = 4>n{s)y n {s), (3) 

leading to a hypergeometric type equation 

a(s)y:(s) + r(s)y' n (s) + Xy n (s) = 0. (4) 
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The first part of the wavefunctions ip n (s) is the solution of the differential equation, 

a(s)<i>'(s) - 7r(s)^(s) = 0, (5) 

where 

T(s)=T(s) + 2n(s), (6) 
and A in (4) is a parameter defined as 

A = A n = -nr'(s) - - {n ~ l) a"(s), n = 0, 1, 2, • • • . (7) 

The t(s) is a polynomial function of the parameter s whose first derivative r'(s) must be 
negative which is the essential condition in choosing the proper solutions. The second part 
of the wavefunctions (3) is a hypergeometric-type function obtained by Rodrigues relation: 

where B n is a constant related to normalization and the weight function p(s) can be found 
by [44] 

a(s)p\s) + [a\s)-r(s)]p( S ) = 0, (9) 
The function n(s) and the parameter A are defined as 



<j'(s)-t(s) I fa'(s) -t(s) x 2 



tt(.s) = ' ' ' ' , ( — r r — J - 3?(a) + fc(7( a ), (f0) 

A = A; + 7r , (s), (ff) 

where n(s) has to be a polynomial of degree at most one. The discriminant under the square 
root sign in Eq. (fO) must be set to zero and then has to be solved for k [44]. Finally, the 
energy eigenvalue equation is simply found by solving Eqs. (7) and (11). 

III. BOUND STATE SOLUTIONS 

The Schrodinger equation for a central molecular potential V(r) can be written as 

(2^ V2 + ~ V{r) ) ^ImM 0) = 0, (12) 

where the representation of the Laplacian operator V 2 , in spherical coordinates, is 
T7 2 -— -— 1 ( 1 9 ■ ft— 1 9 2 \ 
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Here the wave functions ip nlm (r,9,(f)) belong to the energy eigenvalues E ni and V(r) stands 
for the molecular potential in the configuration space and r represents the three-dimensional 
intermolecular distance (X^=i x f) 1 ^ ■ Let us decompose the wave function in (12) as follows: 

ifrraJrM = r~ l u nl {r)Y lm {0^), (14) 

where Yi m (6, (j>) represents contribution from the hyperspherical harmonics that arise in 
higher dimensions. Substituting the wave functions (14) into Eq. (12), the result is [45,46] 

- + | [Era ~ Vir^u^r) = 0, (15) 

where E nt is the bound-state energy of the system under consideration, i.e., E n i < and the 
term is known as the centrifugal term. We also should be careful about the behavior of 
the wave function u n i(r) near r = and r — > oo. Furthermore, u n i(r) should be normalizable 
[47]. 

We can rewrite Eq. (15) for the Hulthen potential as 



d 2 u n i(r) 
dr 2 



+ 



2jjE nl 2/2Ze 2 5 e- &r 1(1 + 1) 



u nl (r) = 0, (16) 



Ti 2 Ti 2 l-e~ Sr r 2 

where E nt is the bound state energy of the system and n and / signify the radial and angular 
quantum numbers, respectively. When I = (s-wave), Eq. (16) with the Hulthen potential 
can be exactly solved [36,48-50], but for the case / ^ 0, Eq. (16) cannot be exactly solved. So 
we must find a new approximation to the entrifugal term to solve the equation analytically. 
The new proposed approximation is based on the expansion of the centrifugal term in a 
series of exponentials depending on the intermolecular distance r and keeping terms up to 
second order. For small 0.4 < 5r < 1.2 [5] (i.e., small screening parameter 5), Eq. (16) 
is very well approximated to centrifugal term. However, for large screening parameter, a 
better approximation to the centrifugal term should be made. Hence, instead of employing 
the usual approximation given in [1-16], we propose an alternative approximation scheme 
casted in the form: 



— w 5 2 [d + v(r) + v 2 (r)] , v(r) = _ ; 



1 

I2 



1 1 

d + -8r—T + 



e 



1 (e Sr - iy 



(17) 
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for the centrifugal term which takes a similar ansatze like the Hulthen potential. Under 
the coordinate transformation r — > x, it is convenient to shift the origin by defining x = 
(r — r )/r , we obtain 

1 1 



d + 



+ 



7 = r 6. 



e 7(l+*) - 1 ' (e7(l+x) _ i) 2 

Further, expanding Eq. (17) around r = r (x = 0), we obtain the following expansion: 



(18) 



1 -2x + 0(x 2 ) = 7 2 (d + 



+ 



7 



+ 



+ 



- 1 ( e 7 - 1) 
2 



- 1 ( e 7 - I) 2 ( e i - 1) 



x + 0(a; 2 ), 



from which we have 



(19) 
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d + 



+ 



- 1 ( e 7 - l) 2 
3 2 



r| — + ^ttp + (^ I=2 ^ (20) 

Therefore, the shifting parameter d is to be found from the solution of the above two 
equations as 

d = - r - — 1 ^ 9 = 0.0823058167837972, (21) 



7' 



- 1 ( e 7 - l)^ 



where e is the base of the natural logarithms, e = 2.718281828459045 and the parameter 
7 = 0.4990429999. 

Therefore, we may cast the centrifugal term as 



lim5 2 

5^0 



1 



-Sr 



+ 



+ 



-Sr 



1 



(22) 



7 2 e7 - 1 (e7 - l) 2 ' 1 - e- dr ' \ 1 - ( 

To conclude, it is important to note that when do = 0, the approximation expression (17) is 
reduced to the usual approximation used in literature [1-16]. The variation of the centrifugal 
term potential 1(1 + l)/r 2 and the proposed approximation expression given in (17) versus 
5r are plotted in Figure 1. Obviously, the approximate centrifugal term potential (17) and 
1(1 + l)/r 2 are similar and coincide in both high-screening as well as in the low-screening 
regimes as shown in Figure 1. 

Inserting the approximation expression (17) into Eq. (16) and changing the variables 
r — > s = e~ 5r through the mapping function s = f(r), where r G [0, oo) or s G [1,0], leads 
us to obtain the following equation 



8(1-8) 



[8(1 - S)Y 



+ (ci - c 2 + 2e^)s - (ci + e 2 nl )s 2 ] u nl (s) = 0, 

(23) 



where 



e nl = xlAE, 



-T272-, AS, = Z(Z + l)do, ci = 2 , c 2 = Z(Z + 1). 



In the present work, we will deal with bound state solutions, i.e., the radial part of the 
wavefunction ip nlm (r,9, 0) must satisfy the boundary condition that u n i(r)/r becomes zero 
when r — > oo, and u n i(r)/r is finite at r = 0. In addition, we require E n \ < ^-AEi, i.e., 
£ni > [36,51-53]. Comparing Eqs. (23) and (2), we obtain the relevant polynomials: 

r( s ) = l-s, a(s) = s(l - s),a(r) = -e 2 nl + (c, - c 2 + 2e?Js - (c, + e 2 nl )s 2 . (25) 

Inserting the polynomials given by Eq. (25) into Eq. (10) gives the polynomial: 

(26) 



7r(s) = -| ± ^V" ,S ' 2 +" ?iN + 



where a = l + 4(ci + £^ — fc), 6 = 4(k — ci + c 2 — 2e^) and c = 4e^. The equation of quadratic 
form under the square root sign of Eq. (26) must be solved by setting the discriminant of 
this quadratic equal to zero: A = b 2 — 4ac = 0. This discriminant gives a new quadratic 
equation can be solved for the constant k to obtain the two roots: 



h,2 = ci - c 2 ± E n i\Jl + 4c 2 . 



(27) 



7r(s) 



(28) 



When the two values of /c given in Eq. (27) are substituted into Eq. (26), the four possible 
forms of 7r(s) are obtained as 

[(e n i - + 4c 2 ) s - e nZ ] for hi = c x - c 2 + V 1 + 4c 2, 
[(e n j + |Vl + 4c 2 ) s - e„i] for fc 2 = ci - c 2 - £„;Vl + 4c 2 . 

One of the four values of the polynomial 7r(s) is just proper to obtain the bound state energy 
states because r(s) given by Eq. (6) has a negative derivative for this value of n(s) [44]. 
Therefore, the most suitable expression of n(s) is chosen as 




7r(s) = 

v ; 2 



1 



:„, + - VI +4c 2 ] s - £ nZ 



(29) 



for A; 2 = c\ — c 2 — e n i VI + 4c 2 . Hence, r(s) and r'(s) are obtained 



r(s) = 1 + 2e ni - 2 



1 



1 + ^+2^1 + 402 



s, 
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r'(s) 



1 



l + e nl + -VI + 4c 2 



(30) 



where r'(s) represents the derivative of r(s). Using Eqs. (25), (29) and (30), the following 
expressions for A and A n are obtained, respectively, 



A = \ n = n 2 + [1 + 2e rd + y/l + 4c 2 ] n, (n = 0, 1, 2, • • • ), 
A = ci - c 2 - \(1 + 2e nl ) [l + VI + 4c 2 J , 



(31) 
(32) 



where n is the number of nodes of the radial wave function u n j(r). When A = A n , an 
expression for e n i is obtained as 

(n + l + 1) 



Cl 



-, (n,/ = 0,l,2,---). 



(33) 



2(n + / + l) 2 

Also, with the aid of Eq. (24), the previous energy equation gives the following bound state 
energy eigenvalue equation: 



jiZe 2 



(n + l + 1) 



-i 2 



h 2 5 (n + l + 1) 



(34) 



In the case of the s-wave (1 = 0), the previous equation turns to be 



h 2 5 2 



E„, = 



fiZe 2 



n+1 



(35) 



2/i \h 2 5(n+l) 2 
which is identical to the ones obtained before using the factorization method [36] , SUSYQM 
approach [3,28,54], quasi-linearization method [55] and NU method [2,6,32]. Further, if we 
take the shift parameter do = in the present approximation, Eq. (34) reduces to 



2 r2 



E n l = 



h z 5 



jj,Ze z 



(n + l + 1) 



(36) 

2/i [h 2 S (n + l + 1) 2 
which is also identical with the energy eigenvalues formula given in Eq. (32) of Ref. [5], Eq. 
(24) of Ref. [43] and Eq. (28) of Ref. [2]. 

Let us turn to the calculations of the wave function y n (s), which is the first part solution 
of hypergeometric-type equation, we need to multiply Eq. (4) by the weight function p(s) 
so that it can be rewritten in self-adjoint form [44] 



[v(s)y' n (s)]' + Xp(s)y n (s) = 0. 



(37) 
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The weight function p(s) that satisfies Eqs. (9) and (37) is found as 



p(s) = s 2£ - l (l-sY 2l+1 \ (38) 

which gives the Rodrigues relation (8): 

Vm{s) = B rd s-^{1 - s)~^-f- [s n+2 ^(l - s) n+2l+1 ] = B nl P^ 2l+l \l - 2s). (39) 

ds n 

Further, inserting the values of a(s) and ir(s) given in Eqs. (25) and (29) into Eq. (5), one 
can find the other part of the wave function as 

(j>(s) = s £ " ; (l-s) (m) . (40) 

Hence, the wave functions in Eq. (3) become 

u nl {s)=M nl s^{\ - s) l+1 P^ 2l+l \l - 2s), s e [1,0) (41) 

where M n i is the normalization constant to be determined in the next section. Finally, the 
unnormalized radial wave functions are obtained as 

Vw(r,0, 0) = Mmr- 1 (e~ 5r ) £nl (l- e - Sr ) l+1 F 1 (-n, n+2 (e nl + 1 + 1); 2e nl + l; e- 5r )Y lm (0, 0). 

(42) 

Thus, the Jacobi polynomials can be expressed in terms of the hypergeometric functions 
[56]: 

Pi a ' b \l -2x) = 2 Fi(-n,ra + a + 6 + l;a + l;a;), (43) 

oo 

where 2 Fi(a, b; c; x) = r ^\ b) ^ ^tillk)^ i~ Tne hypergeometric function 2 F 1 (a, b; c; x) 

k=0 

is a special case of the generalized hypergeometric function [56] 

p F q (a u a 2 , ■ • • , a p ; f3 1: • • • , f3 q ; x) = ^ >k Jk ^(-, (44) 

fc=0 {Pl)k \P2)k ■ ■ ■ \Pq) K - 

where the Pochhammer symbol is defined by (y) k = T(y + k)/Y(y). 
In the case / = 0, the above wave functions become 

^ n (r) = D n r~ l (e- 5r f" (1 - e'^F^-n, n + 2(e n + l); 2e n + 1; e' 5r ), (45) 

with e n = ffa^n+i) ~ r± ~2~ an< ^ D n is another normalization factor. This result is consistent 
with the NU method [2]. Further, if we take the shift parameter d = in the present 

10 



approximation, Eq. (42) reduces to the form 

i> nlm (r,e, 0) = D nl r~ l (e~ 5r ) £nl (l-e' Sr ) l 2 +1 F 1 (-n, n+2 (e nl + 1 + 1); 2e nl + l; e- 5r )Y lm (0, 0), 

(46) 

with e n i = tfs(n+i+i) ~ n+l 2 +1 an< ^ ^ni is a normalization factor. The critical screen- 
ing = tff^ptf at which E ni = has wave functions: ijj nlm (r,6,(j)) = D n ir~ l {l - 
e s cr y + ip(o,2i+i)^ _ 2e -^)Y lm {e,(t>)- 

In order to show the accuracy of our analytical results, we present the numerical data in 
support of the results obtained in Eqs. (34) and (42) which are the main analytic results 
obtained in this work. Therefore, we calculate the energy eigenvalues for Z — 1, n and I ar- 
bitrary quantum numbers as a function of the screening parameters S. The results calculated 
in Tables 1 and 2 by using Eq. (34) are compared with those obtained with the help of the 
numerical integration [37], asymptotic iteration [5], variational [37], SUSY [3] and the re- 
cently proposed approximation [43] methods. Tables 1 and 2 show that our results obtained 
with the new approximation scheme with the NU method are in high agreement with those 
obtained by numerical integration method [37] for short potential range (small screening 
parameter 5). However, the slight differences in the energy eigenvalues from the numerical 
integration method [37] are observed for long potential range (large screening parameter 
5). Therefore, our approximated numerical results are closer to the numerical integration 
results than the results obtained via AIM [5] using Eq. (32) and also the recently proposed 
approximation scheme [43] using Eq. (19) for small and large screening parameter S values. 
Thus, the present approximation form (22) to the centrifugal term can highly improve the 
accuracy of calculating the energy eigenvalues for the Hulthen potential than the recently 
proposed approximation (5) given in [43]. 

IV. NORMALIZATION OF THE RADIAL WAVE FUNCTION 

Using s(r) = e~ Sr and Eq. (41), we are able to express normalization condition 

J °° u{r) 2 dr = 1 as 

f 1 s 2e nl -l {1 _ s) 2l + 2 [p^ nl ,2l + l) {1 _ 2s) ^ ds=h (47) 
° JO 

Unfortunately, there is no formula available to calculate this key integration. Neveretheless, 
we can find the explicit normalization constant N n i. For this purpose, it is not difficult to 
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obtain the results of the above integral by using the following formulas [56-59] , 

1 fx — 1 \ f X + 1^ 



P^\x) = (n + a)\(n + (3)\J2 



o p\(n + a -p)\ ((3 + p)\ (n + p)\ \ 2 J \ 2 



and 



B(x,y) = - ty-'dt = , Re(x),Re(y) > 0. 

Thus, the normalization constant M n i is now obtained as 



(48) 
(49) 



nl 



{n + 2l + l)\T(2e nl + n + l) 



i 



ST(2e n i + 2n + 2l + 4) 

n ■ 

T{2e nl + 2n + l) E UvMp*) 

p,q=0 



where 



U = (-l) p p\r(2e nl +n-p + l)(2l+p + l)\(n + p)\, 
f q = (-iyq\T(2e nl + n-q + l)(2l + q + l)\(n + q)\, 

f Ptq = (2l + p + q + 2)\. (50) 

V. CONCLUSIONS 

In this work, we have proposed an alternative improved approximation scheme for the 
centrifugal term and used this approximation scheme together with the NU method to solve 
the Schrodinger equation with any orbital angular momentum number I for the Hulthen 
potential. The bound state energy eigenvalues and the unnormalized radial wavefunctions 
have been calculated in analytical and numerical way. The analytic expressions for the energy 
eigenvalues and wavefunctions have been reduced to the s-wave case and the do = case 
(usual approximation) [1-16]. Our numerical results obtained by the approximation scheme 
given in expression (22) for the centrifugal term has been found to be more effective than the 
numerical results of the recently proposed approximation (5) of Ref. [43] and the commonly 
used approximation in generating the energy spectrum of the Hulthen potential. Our results 
in Tables 1 and 2 for small screening 5 values show that the present approximation is in high 
agreement with the numerical integration and variational methods [37] whereas it is in quite 
good agreement for large screening 5 values. The present approximation method is simple, 
practical and powerful than the other known methods [2,5,43]. This new method can be used 
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for many quantum models to improve the accuracy of energy eigenvalues for few potential 
models of the exponential-type like the hyperbolical and Manning- Rosen potentials (cf. e.g., 
Refs. [60,61].) 
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FIG. 1: A plot of the variation of the centrifugal term, 1/r 2 and its corresponding approximation 



expression S 2 



versus 5r, where the screening parameter 5 changes from 5 = 0.050 



(e^l) 2 _ 

to S = 0.250 in steps of 0.050. The parameters are in atomic units (h = = e = 1) and Z = 1. 



17 



TABLE I: Bound energy spectra of the Hulthen potential as a function of the screening parameter 
5 for 2p, 3p and 3d states for Z = 1 in atomic units (h 2 = jj, = e = 1) 

State 5 Present" Previous [43] b Numerical [37] AIM [5] Variational [37] SUSY [3] 



0. 


.025 


0. 


.1127611 


0.1126344 


0.1127605 


0. 


1128125 


0. 


.1127605 


0. 


.1127605 


0. 


.050 





.1010442 


0.1009128 


0.1010425 


0. 


1012500 


0. 


.1010425 


0. 


1010425 


0. 


.075 


0. 


.0898495 


0.0898350 


0.0898478 


0. 


.0903125 


0. 


.0898478 


0. 


.0898478 


0. 


100 





.0791769 


0.0794011 


0.0791794 


0. 


.0800000 


0. 


.0791794 


0. 


.0791794 


0. 


.150 


0. 


.0593981 


0.0604650 


0.0594415 


0. 


.0612500 


0. 


.0594415 


0. 


.0594415 


0. 


.200 


0. 


.0417078 


0.0441045 


0.0418860 


0. 


.0450000 


0. 


.0418860 


0. 


.0418854 


0. 


.250 


0. 


.0261059 


0.0303195 


0.0266111 


0. 


0312500 


0. 


.0266108 


0. 


.0266060 


0. 


300 


0. 


.0125925 


0.0191101 


0.0137900 


0. 


.0200000 


0. 


.0137878 


0. 


0137596 


0. 


350 


0. 


.0011675 


0.0104763 


0.0037931 


0. 


.0112500 


0. 


.0037734 


0. 


.0036146 


0. 


.025 


0. 


.0437072 


0.0436848 


0.0437069 


0. 


.0437590 


0. 


.0437069 


0. 


.0437068 


0. 


.050 





.0331623 


0.0332390 


0.0331645 


0. 


.0333681 


0. 


.0331645 


0. 


0331632 


0. 


.075 


0. 


.0239207 


0.0242183 


0.0239397 


0. 


.0243837 


0. 


.0239397 


0. 


.0239331 


0. 


.100 





,0159825 


0.0166227 


0.0160537 


0. 


.0168056 


0. 


.0160537 


0. 


.0160326 


0. 


.150 


0. 


.0040162 


0.0057067 


0.0044663 


0. 


.0058681 


0. 


.0044660 


0. 


.0043599 


0. 


.025 


0. 


.0436044 


0.0435371 


0.0436030 


0. 


.0437587 


0. 


.0436030 


0. 


.0436030 


0. 


.050 





.0327508 


0.0329817 


0.0327532 


0. 


.0333681 


0. 


0327532 


0. 


0327532 


0. 


.075 


0. 


.0229948 


0.0238893 


0.0230307 


0. 


.0243837 


0. 


.0230307 


0. 


.0230306 


0. 


100 





.0143364 


0.0162600 


0.0144842 


0. 


.0168055 


0. 


.0144842 


0. 


.0144832 


0. 


150 


0. 


.0003124 


0.0053907 


0.0013966 


0. 


.0058681 


0. 


.0013894 


0. 


.0132820 



°The present approximation. 
b Jia et al approximation. 
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TABLE II: Bound energy spectra of the Hulthen potential as a function of the screening parameter 
8 for 4p, 4d, 4f, 5p, 5d, 5/, 5g, 6p, 6d, 6/ and 6g states for Z = 1 in atomic units (Ti 2 = n = e = 1) 
State 5 Present" Previous [43] b Numerical [37] AIM [5] Variational [37] SUSY [3] 



4p 


0. 


.025 


0, 


.0199486 


0.0199625 


0.0199489 


0. 


0200000 0.0199489 


0. 


.0199480 




0. 


.050 


0. 


.0110442 


0.0111938 


0.0110582 


0. 


0112500 0.0110582 


0. 


0110430 




0. 


.075 


0, 


.0045370 


0.0049439 


0.0046219 


0. 


.0050000 0.0046219 


0. 


.0045385 




0. 


.100 





.0004269 


0.0012128 


0.0007550 


0. 


.0012500 0.0007532 


0. 


.0004434 


4d 


0. 


.025 


0. 


.0198457 


0.0198877 


0.0198462 


0. 


0200000 0.0198462 


0. 


.0198460 




0. 


.050 


0. 


.0106327 


0.0110819 


0.0106674 


0. 


0112500 0.0106674 


0. 


.0106609 







.075 


0. 


.0036111 


0.0048327 


0.0038345 


0. 


.0050000 0.0038344 


0. 


.0037916 


4/ 


0. 


.025 


0. 


.0196914 


0.0197756 


0.0196911 


0. 


.0200000 0.0196911 


0. 


.0196911 







.050 





.0100154 


0.0109150 


0.0100620 


0. 


.0112500 0.0100620 


0. 


.0100618 







.075 


0, 


.0022222 


0.0046682 


0.0025563 


0. 


.0050000 0.0025557 


0. 


.0025468 


hp 


0. 


.025 





.0094017 


0.0094325 


0.0094036 


0. 


.0094531 


0. 


0094011 




0. 


.050 


0. 


.0026067 


0.0027900 


0.0026490 


0. 


.0028125 


0. 


.0026056 


5d 





,025 





,0092988 


0.0093914 


0.0093037 


0. 


.0094531 


0. 


.0092977 




0. 


,050 


0. 


,0021952 


0.0027454 


0.0023131 


0. 


.0028125 


0. 


.0022044 


5/ 





.025 





.0091445 


0.0093898 


0.0091521 


0. 


.0094531 


0. 


.0091507 







.050 


0, 


.0015779 


0.0026791 


0.0017835 


0. 


.0028125 


0. 


.0017421 


5g 


0. 


.025 


0, 


.0089387 


0.0092480 


0.0089465 


0. 


.0094531 


0. 


.0089465 







.050 





.0007549 


0.0025920 


0.0010159 


0. 


.0028125 


0. 


.0010664 


Qp 


0. 


.025 


0. 


.0041500 


0.0041899 


0.0041548 


0. 


.0042014 





.0041493 


6d 





,025 





,0040471 


0.0041671 


0.0040606 


0. 


.0042014 


0. 


.0040452 


6/ 


0. 


,025 


0. 


.0038927 


0.0042014 


0.0039168 


0. 


.0042014 


0. 


.0038901 


6g 





.025 





.0036870 


0.0040876 


0.0037201 


0. 


.0042014 


0. 


.0036943 



a The present approximation. 
6 Jia et al approximation. 
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